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Abstract. Within the framework of Bonnor's exact solution describing a massive 
I magnetic dipole, we study the motion of neutral and electrically charged test particles. 

^Jpj In dependence on the Bonnor spacetime parameters, we determine regions enabling 

i— i the existence of stable circular orbits confined to the equatorial plane and of those 

levitating above the equatorial plane. Constructing Poincare surfaces of section and 
j>. recurrence plots, we also investigate the dynamics of particles moving along general 

t-H off-equatorial trajectories bound in effective potential wells forming around the stable 

circular orbits. We demonstrate that the motion in the Bonnor spacetime is not 
integrable. This extends previous investigations of generalized Stormer's problem into 
the realm of exact solutions of Einstein-Maxwell equations, where the gravitational 
and electromagnetic effects play a comparable role on the particle motion. 
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1. Introduction 

Since the seminal work by Stormer [TJ, many authors have investigated properties of 
motion of electrically charged particles through a dipole-type magnetic field as a model 
for magnetospheres surrounding planets. The original motivation of the early works was 
to understand the trajectories and radiation of plasma particles which arrive with the 
solar wind and become captured by the Earth magnetic field. Indeed, radiation belts are 
now known to be formed by ions and electrons whose motion is governed by magnetic 
forces. Beside these classical studies, the Stormer's work was used in the context of 
General Relativity to study the stable orbit around Kerr black hole by Ruffini and 
Wheeler. Their work is reproduced among others in figure 18 of the book [2] and has 
been quoted in [3] in paragraph 104. All the definitions of the stable and unstable orbits 
were there basically obtained by using the Stormer's approach in General Relativity (see 
also (4|). 

Later on, different generalizations of the original Stormer's discussion have been 
developed in order to describe various aspects of the problem and, to capture more 
complicated situations (5j[6]. For example, in a case of dust particles, only a weak 
electric charge is established on the surface of the particles, i.e., small charge to mass 
ratio develops by photoionization or selective charge accretion [7] . Therefore, one needs 
to take also gravitational force into account. At the same time the influence of rotation 
induced electric field becomes non-negligible and has to be combined with the magnetic 
component. This is the essence of the so-called generalized Stormer problem [5j|8]. 
Total gravitational and electromagnetic forces then determine regions of stable motion, 
in particular, whether off-equatorial potential wells (halo lobes) develop where particles 
can be captured in permanent circulation around the central body. The problem is 
relevant in the context of accreting compact objects, where the motion within the halo 
lobes (halo motion) can describe the overall global motion of plasma particles through 
corona of an accretion disc. 

Relativistic effects change the regions of stability near a compact star or a black 
hole. In order to reveal the existence of the halo motion under strong gravity regime, we 
adopt the formalism of General Relativity. To this end, a slowly rotating compact star 
endowed with an aligned magnetic dipolar field, or a magnetized black hole embedded 
in an external magnetic field of a current loop outside the horizon were previously 
studied (9 11 . Furthermore, an asymptotically uniform Wald magnetic field parallel 



to the rotation axis of Kerr black hole was also considered 10 12 . Regions of stable 
circular halo motion were found for suitable combinations of model parameters, such 
as the generalized energy and angular momentum of the particles, their specific electric 
charge, the magnetic dipole moment of the central body (defining the intensity of the 
large scale magnetic field), and angular velocity of rotation. Note that an approximative 
analysis of the halo motion in strong gravitational fields can be also formulated by 



employing the pseudo- Newtonian approach (see, e.g., paper 13 for an overview and 



references). Although this method does not allow a mathematically rigorous discussion 
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of the problem, it can provide useful clues about the effects and, indeed, it was found [9] 
that the existence and position of the stable halo circular orbits of charged particles are 
reproduced very accurately by the Paczyhski & Wiita model. 

Whilst the above-mentioned examples include gravitational effects by adopting the 
metric of a black hole, the electromagnetic influence has been so far considered in terms 
of only a weak (test) field approximation. The Kerr-Newman black hole offers a way 
to study the problem within the framework of an exact solution of combined Einstein- 
Maxwell equations. However, this particular example turns out to be quite uninteresting, 
because the halo lobes were not found above the outer horizon [9]. The lack of stable 
circular halo orbits outside the horizon of Kerr-Newman black hole is not very surprising 
(it can be ascribed to the special character of this solution), nonetheless, we want to 
explore the occurrence of these orbits and study the related halo motion in an exact 
solution that can be used to describe a massive magnetic dipole as the closest analogy 
to the Newtonian generalized Stormer problem. To do this, we employ the Bonnor 
spacetime as the suitable example of such a background. 

Test particles motion in strong gravitational and electromagnetic fields has been 
ranked among basic problems in relativistic astrophysics, widely studied since the 
conception of General Relativity, still having high theoretical, practical and educational 



significance (see, e.g., the recent nice serie of papers 14 16 or the paper 17 and the 
reference therein). 

2. Bonnor's solution and equations of motion 

Bonnor derived a special solution 08] of combined Einstein-Maxwell equations 
describing a static massive source carrying a magnetic dipol. In (t, r, 9, <fi) coordinates 
and c = G = 1 units the solution reads 

,2 f p Y , 2 P 2 Y\, 2 ^ 1/l2 , Y 2 Z sin 2 9 lj2 
ds 2 = - (^J dt 2 + Q^(dr 2 + Zd9 2 ) + p2 d0 2 , (1) 

ur sin 2 9 . 

A a =(0,0,0 / - ), (2) 
where we adopted the notation 

P = r 2 - 2ar - b 2 cos 2 9, (3) 

Q = ( r - a) 2 - (a 2 + b 2 ) cos 2 9, (4) 

Y = r 2 - b 2 cos 2 9, (5) 

Z = r 2 - 2ar - b 2 . (6) 

The solution is asymptotically flat (for a = exactly flat) and of Weyl tensor type 
[1111]. It is characterized by two independent parameters a and b, related to the total 
mass of the source M = 2a (from the g tt asymptotic behaviour) and to the magnetic 
dipole moment \x = 2ab (from the asymptotic behaviour). The Bonnor line element 
((!]) manifests relatively complicated singular behviour at P = 0, Q = 0, Z = and 



Y = 0, nature of which is described in details in the paper 19 . However, here we are 
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interested in the regular part of the spacetime Z > only, limited by the horizon at 
r h = a + \/ a 2 + b 2 . 

The Bonnor field can be interpreted as a superposition of two magnetic poles 
of strengths —a and a, both with mass a, that are situated on the symmetry axis 
and separated by the coordinate distance 2b [l9|jf] The problem how to keep this 
magnetic monopoles appart in an equilibrium was answered, e.g., by adding an external 
magnetic field into the solution [20). Such a configuration is analogous to two oppositely 
electrically charged Reissner-Nordstrom masses, which can survive in equilibrium due 
to an induced electric field mechanism, as recently confirmed within the framework 



of an exact solution of the static Einstein-Maxwell equations 21,22]. A sophisticated 



analysis of the electrostatic two-body problem has been recently presented also in 23 24 



completed by the investigation of electric field lines, exhibiting the Meisssner effect found 



in such system 25 



Naturally, a number of more recent exact solutions have been discovered, which can 



represent a massive magnetic dipole in relativity (see, e.g., papers 26-29]). They differ 
from each other by the structure of the field near the central body, where singularities 
can occur. But despite a large number of them, Bonnor's solution Il8] plays a prominent 
role thanks to its simplicity and explicit form that allows us to interpret the solution 
and describe its properties in great details. The metric Q actually represents a 



magnetostatic limit of a more general exact solutions 30,31] suggested to describe the 
exterior field of a rotating neutron star. In the case of Bonnor's solution, the rotation 
is not considered and the value of a quadrupole mass moment is fixed by values of the 
mass and dipole parameters a and b 

Q = ^--M* = 2a{b*-a% (7) 

vanishing and changing its sign for \b\ = a. These features support an exclusion of 
the pure Bonnor spacetime from astrophysically realistic ones (which are not static, 
instead they are dragged into a rotation by spinning source, often also keeping the 
quadrupole moment as another independent parameter). Thus, the investigation of 
particle motion in the Bonnor's static magnetic dipole background presented here has 
mainly a theoretical importance; for example, it allows us to study the onset of chaos in 
General Relativity. From the theoretical point of view, however, it is interesting to point 
out that the frame-dragging, being typically related with spacetimes of spinning sources, 
can be also found in the vicinity of static source endowed with an electromagnetic 
structure 32 , 33 . In more details, as shown in the paper 33 studying a rotating 
magnetic dipole bearing an electric charge, the vorticity tensor is nonvanishing and 
the frame-dragging occurs even in the special case when the angular momentum of the 
source is zero. This is due to the presence of the magnetic dipole moment and the 

| In the case of b = the proper distance between the magnetic poles is nonzero and the geometry 
does not tend to the Schwarzschild one. 
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additional electric charge of the source, which is, however, not the case of the Bonnor's 
massive dipole studied here. 

According to the standard approach to investigation of particle dynamics, we 
construct the super-Hamiltonian [4] 

n = Ig^fa - qA^ v - qA„), (8) 

where q is the charge of particle, 7r M is the generalized (canonical) momentum, is the 
metric tensor (TXT) and denotes the electromagnetic vector potential ([2]). The latter 
is related to the electromagnetic tensor by F^ u = A v>li — A^ v . The related Hamiltonian 
equations which we integrate to obtain general trajectories of particles are given as 



dx 11 



p F 



&H dir. 



dX 



on 



(9) 



where A = r/m is the affine parameter, r denotes the proper time, m stands for the 
rest mass of the particle and p 11 is the standard kinematical four-momentum for which 



the first equation reads p u 



qA>*. 



The Bonnor's solution obeys two obvious symmetries defined by Killing vectors 
fields if" = 5i and £ M = 5^, which ensure the existence of two related constants of 
motion, quaranteed also by second Hamiltonian equation. These are the specific (per 
unit mass) energy and the specific angular momentum 

2 



E 



u t - qA t 



L — u ( j, + qA^ = Z sin 2 9 



dt 
dr 
Y 



(10) 



dr 



[ID 



where q = q/m is the specific charge of particle. 

The motion in the meridian plane can be described very well by the equation 



E 



P 4 
Q 3 Z 



( dr 
[dr 



with the effective potential given by the relation 



V 2 



p2 

Y 2 



1 + 



P 2 



Y 2 Z sin 2 6 



L 



qH- 



r sin 6 



P 



(12) 



(13) 



Because of the properties of the metric (squared b parameter) and the effective potential 
(b parameter in product with q, squared L—qA^ term), the effective potential is invariant 
under combinations of simultaneous sign reversals maintaining sgn(Lg), sgn(L6) or 
sgn(g6). Since the metric is static and axially symmetric, we expect that closed, stable 
planar orbits should be possible in the equatorial plane (9 = 7r/2). The self-consistent 
magnetic field obeys the same symmetries in this solution, thus we can also expect the 
possibility of the off-equatorial azimuthal circulation of electrically charged matter when 
the Lorentz force is in equilibrium with gravitational and centrifugal forces. 
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The circular orbits of constant radius and latitude (circular halo orbits) are 
determined by stationary points of the effective potential, i.e., by the conditions 

d r V eS = 0, d e V eS = 0. (14) 

Being interested in the stability of the motion against the perturbations, we have to 
construct the Hessian matrix 

H _ I d ?Ve& d r d e V eS \ 

\d e d r v cS d*v eS )■ [Lt>) 

Circular orbits at r = r c and 9 = 9 C stable against all the perturbations (minima of V e s) 
must satisfy the conditions 

d^V eS \ r=rc; = c > 0, detH\ r=rct e = e c > 0. (16) 

In particular, the stability of circular orbits against radial perturbations (important 



mainly in the case of equatorial motion) is guaranteed by the first from conditions (16). 
On the other hand, fully unstable circular orbits (maxima of V e g) satisfy the conditions 

d*V eS \ r=rct e =0c < 0, det#| r=rc) =(?c > 0. (17) 



Investigating the circular motion, it is also useful to monitor the azimuthal component 
of the linear velocity 

= y/=gT t (L - q Atj,) 



v <t> = ^i = B 7= • 



3. Circular equatorial motion 

The equatorial circular motion of charged particle at radius r is characterized by the 
specific angular momenta 

qab{2F + rZ) ± rZ y / (2aF + (qab) 2 ) 



(r - 2a)F 



(19) 



following from the first of conditions (14) or, alternatively, by the related azimuthal 
component of the linear velocity 

v c ± — v(L — L c± , E — V eS ). (20) 

In the special case of uncharged (q = 0) particles, the reality condition F = 



Aab +r(r — 5a) (r — 2a) > in relation (19) restricts the existence of equatorial circular 



orbits by the photon circular orbits at radii given by the solutions of F = 0, in the 
(r x 6) plane determined by the curve defined by the functions 



r(5a - r)(r - 2a) 

V = ±V Ta • (21) 

Extrema of these functions are located at r = 3.786 a, taking values b = ±1.433 a. In 
the case of b = 0, the inner photon orbit coalesces with the horizon at r^ = 2 a, the 



outer one is located at r = 5 a. As it follows from the investigation of conditions (16) 
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with 6 C = 7r/2, in the (r x b) plane, the r-region of fully stable orbits and orbits unstable 
against radial perturbations are separated by the curve defined by the functions 

^ , (±r(r - 2a) v/(3r 2 - lOar + 6a 2 ) (3r 2 + 22ar - 26a 2 ) 

~ _± l 32^ (22) 



II1S 



3r(r 3 - 18a 2 r + 28a 3 



1/2 



32a 2 

maximum and minimum of which are located at r = 7.236 a and take values b = 
±2.115 a. The curves and b p ^ interstect each other just in the extrema of 6 p h- Fully 
stable orbits and orbits stable against radial perturbations (but not fully stable) are 
separated by the curve defined as 



«- = ±r ^ra- (23) 

We present stability regions and the introduced functions in figure [T] along with 
the typical behaviour of the specific angular momentum L° c+ = L c+ (q = 0) and velocity 
v® + = v c+ (q = 0), plotted for different values of b from intervals defined by the maxima 
of 6 p h and 6^ s . Note that positions of minima of L° c+ separate the r-regions of stability; 
the descending parts correspond to r-stable orbits and ascending parts to r-unstable 
orbits. For b = 0, we find the stable orbits limited from below by r ms = 11.359 a. 
The specific angular momentum L° c+ diverges at the radii of photon orbits (related 
velocity reaches the speed of light here) and vanish at the radius of horizon as 
well as the velocity. Because of the absence of the electric charge, we have here two 
specific angular momentum profiles L° + and L° c _ (19), differing only in the sign, and the 
effective potential resistant to a change of the sign of L. Consenquently, the stability 
regions illustrated in figure [TJ calculated for L = L° c+ profile, remains the same for 
L°_ profile, only the signs of specific angular momentum and velocity contours will be 
different. 



In the general case of charged particles g / 0, we get the functions (22) and 
(23) heavily dependent on the specific charge, i.e., we have 6 ms = 6 ms (r;a, q) and 
6 mcs = 6 mcs (r; a, q). The relevant expressions, however, are rather long to be mentioned 
here and the analytical discusion is very complex. Nevertheless, in figure [2j we plot the 
representative behaviour of these functions for the fixed value q = 4.716, giving us one of 
the most complex configurations of the stability regions. The stability regions shown in 
figure [2] are constructed providing L = L c _ and, will change when L = L c+ . But because 
of the symmetries of the effective potential, just in the way of inverting along the 6 = 
axis and changing the signs of specific angular momentum and velocity contours. In 
figure |4| we plot the stability regions from this class as well, i.e., corresponding to the 
L = L c+ root, and setting q = 8. 

Positions of minima of L c ± separate the r-regions of stability, the descending parts 
correspond to r-stable orbits and ascending parts to r-unstable orbits, as well as in the 
uncharged case. The situation here, however, is generally much more interesting. At 
first, note that thanks to the electric charge, we can find also fully unstable circular 
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Figure 1. Left: Regions of stability for circular orbits of uncharged particles moving 
in the equatorial plane with the specific angular momentum profile given by L® + : region 
of orbits stable against all the perturbations (white), region of orbits stable against 
the radial perturbations but not fully stable (yellow R-I) limited by the curve b^ cs and 
region of orbits unstable against the radial perturbations (blue R-II) limited by the 
curve 6^ s . The circular orbits do not exist in the oval-like black region limited by the 
curve b p h denoting positions of circular photon orbits; the triangle-like black regions 
represent the regions under the horizon Z = 0. The dashed (k = 5, 1 = 7, m = 9) and 
solid (K = 0.4, L = 0.6) curves denote contours of specific angular momentum L° + (in 
units of a) and related velocity The horizontal dotted lines and their intersections 
with the stability regions and limiting curves suggest the behaviour of L® + and v® + as 
shown in the right panels of this figure. Right: Specific angular momentum L®, and 
related velocity profiles for 6 = (solid), b = 1.3a (dot-dashed), b = 1.7 a (dashed) 
and b = 2.5 a (dotted). The angular momentum diverges at the radii of photon orbits, 
where the related velocity reaches the speed of light, and vanishes at the radius of 
horizon as well as the velocity. 



orbits, satisfying conditions (17). Moreover, we have two completely different specific 
angular momentum profiles L c+ and L c _ (19) here. Thus, for instance, for a special 
combination of the spacetime parameter b and specific charge q, the L c+ profile diverges 
at the radii of photon orbits (related v c+ reaches the speed of light here) and vanishes 
at the radius of horizon as well as v c+ . We call the related circular orbits Larmor orbits. 
On the other hand, the L c _ profile takes finite values at the radii of photon orbits where 
the related w c _ are sublight. The related orbits are then anti-Larmor. Apparently, the 
circular orbits of charged particles can exist in the region between the photon orbits as 
well. Such a situation is well illustrated in figure [5j 
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Figure 2. Regions of stability in two different zooms for circular orbits of charged 
particles with q — 4.718 moving in the equatorial plane with the specific angular 
momentum profile given by L c _: region of orbits stable against all the perturbations 
(white), region of orbits stable against the radial perturbations but not fully stable 
(yellow R-I) limited by the curve 6 mcs , region of orbits unstable against the radial 
perturbations (blue R-II) limited by the curve o ms and region of orbits unstable against 
all the perturbations (red R-III) . The circular orbits do not exist in the oval- like black 
region; the triangle-like black regions represent the regions under the horizon Z = 0. 
The dashed (k = -12, 1 = -4, 1* = 4, k* = 12) and solid (K = -0.4, L = -0.6, M 
= —0.995, M* = 0.995) curves denote contours of specific angular momentum i c _ (in 
units of a) and related linear velocity t> c _. The thick solid oval-like curve 6 p h denotes 
positions of circular photon orbits. The thick dashed curve denotes the zero specific 
angular momentum contour. The horizontal dotted line b — —3a and its intersections 
with the stability regions and limiting curves suggest the behaviour of L c _ and v c _ as 
shown in figure [3] 




10 15 20 25 30 35 
r[a] 



Figure 3. Specific angular momenta L c _ (solid) and L c+ (dot-dashed) and related 
linear velocities v c - (solid) and u c + (dot-dashed) profiles for particles with specific 
charge q — 4.718 moving in Bonnor spacetime with b = —3 a. The angular momenta 
profiles coalesce with each other at the radius of horizon, where the related velocities 
vanish. 




Figure 4. Regions of stability in two different zooms for circular orbits of charged 
particles with q — 8 moving in the equatorial plane with the specific angular momentum 
profile given by L c+ : region of orbits stable against all the perturbations (white) , region 
of orbits stable against the radial perturbations but not fully stable (yellow R-I) limited 
by the curve 6 mcs , region of orbits unstable against the radial perturbations (blue R-II) 
limited by the curve 6 ms and region of orbits unstable against all the perturbations (red 
R-III). The circular orbits do not exist in the oval- like black region; the triangle-like 
black regions represent the regions under the horizon Z = 0. The dashed (k= —12, 1 
= -4, 1* = 4, k* = 12) and solid (K = 0.4, L = 0.6, M = 0.995, M* = -0.995) curves 
denote contours of specific angular momentum L c+ (in units of a) and related linear 
velocity w c+ . The thick solid oval- like curve & p h denotes positions of circular photon 
orbits. The thick dashed curve denotes the zero specific angular momentum contour. 
The horizontal dotted line b — 0.42 a and its intersections with the stability regions 
and limiting curves suggest the behaviour of L c+ and v c+ as shown in figure [5] 



4. Circular halo motion 



Loci of circular halo orbits (9 ^ vr/2) determined by conditions (14) are characterized 
by the specific angular momenta and charges 

— Y sin 6 (r 2 + b 2 )x/ra , . 

A/ZA(r 2 + b 2 cos 2 6) 



aY 2 S 2 

~2ab ]J rZA sin 2 9{r 2 + b 2 cos 2 
where 

S = r 4 - 2ar 3 - 3b 2 r 2 sin 2 9 + 2ab 2 r{2 sin 2 9 - 1) - b 4 cos 2 9, (26) 



A = 3r 4 - 10ar 3 - 26 2 r 2 cos 2 9 - 6ab 2 r cos 2 9 - b 4 cos 4 9. 



(27) 
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Figure 5. Specific angular momenta L c+ (solid) and L c _ (dot-dashed) and related 
linear velocities v c+ (solid) and u c _ (dot-dashed) profiles for particles with specific 
charge q — 8 moving in Bonnor spacetime with b = 0.42 a, plotted in different zooms. 
The L c+ profile diverges at the radii of photon orbits (vertical solid lines), where the 
related velocity reach the speed of light, and coalesce with L c _, among others, at the 
radius of horizon, where the related velocities vanish. The L c _ profile takes finite 
values at the radii of photon orbits, where the related velocity is sublight. 



The related linear velocity can be specified as 

v h ± = v(L = L h ±, q = q h ±, E = V eS ). (28) 

Apparently, the existence of circular halo orbits in the Bonnor spacetime in the region 
Z > is determined by the condition A > 0. Next, we are mainly interested in the 
stable motion, i.e., in the orbits which satisfy conditions (16) for the local minima of 
the potential. The regions of existence and stability of the circular halo orbits are 
illustrated for the representative latitude 6 = it/ 3 in figure [6j providing that L = 
and q = It is important to point out that the regions of existence and stability 
remain the same after taking the other roots and g^- into account for conditions 
(16) analysis. Moreover, the figure appears qualitatively the same after changing 9. We 
also plot behaviour of the angular momenta, related velocity and specific charges profiles 
in figure [7j 

Apparently, the special geometry of the Bonnor spacetime enables the existence 
of the circular halo orbits of uncharged particles (circular halo geodesies) as well. 
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2 4 6 8 



r[a] 

Figure 6. Regions of stability for circular halo orbits of charged particles with the 
specific angular momenta and charges profiles given by Lh+ and <?h+ : region of orbits 
stable against all the perturbations (white), region of orbits with an instability (blue 
R-I) and region with no circular halo orbits (gray R-II). The black regions denote the 
area under the horizon Z = 0. The solid bow-like (more vertical) curves (K = —0.7, L 
= —0.5, L* = 0.5, K* = 0.7) denote contours of linear velocity 1^+, the dashed curves 
(k = —6, 1 = —2.5, k* = 6) denote the contours of specific angular momentum ih+ (in 
units of a) and the solid (more horizontal) curves (a = —0.7, f3 = — 7 a* = 0.7, (3* = 7) 
denote the contours of specific charge qh+ ■ The thick solid curve separating the regions 
RT and R-II corresponds to part of the A = curve; the transversal thick red curves 
correspond to the S = curve, denoting positions of circular halo geodesies (dividing 
the positive and negative contours of Lh+ and t>h+); the positions of circular halo 
photon orbits correspond to their intersections (red bullets). The horizontal dotted 
lines b = 3.5 a, b = 2 a and b = 1.668 a denote the values for which the specific angular 
momenta, linear velocities and specific charges profiles are illustrated in figure [7] 



As it follows from relations d r V e s\ q=0 = and dgV e s\ q= o = after eliminating L or, 
directly from relation (25), their positions are determined by the condition S = 0. The 
angular momentum of the uncharged particles moving along circular halo geodesies are 
determined by the relation L° ± = Lh±/sgn(S') (where the radii and latitudes are related 
by the formula 5 = 0). The relation for L^± diverges for A = 0, therefore the condition 
A = together with S — enable us to determine positions of even circular halo photon 
orbits, as illustrated in figure |HJ The reality conditions here restrict their existence to 
the regions r > 3.333 and for spacetimes with |fe| > 1.361. 
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Figure 7. Specific angular momenta Lh+ (solid) and Lh- (dot-dashed), specific 
charges q\ 1+ (solid) and qh- an d related linear velocities «h+ (solid) and i>h- for particles 
moving in Bonnor spacetimes with b = 3.5a, 6 = 2a and b = 1.668a. The angular 
momenta and charge profiles diverge at the radii of horizon (vertical solid line) given 
by Z = or at the radii given by A = (vertical dashed line) , if they are above the 
horizon. Of course, at the radii given by S = (determining uncharged circular halo 
orbits) the charges vanishes. The velocities reach the speed of light at the radii given 
by A = 0. 



5. Regular and chaotic dynamics 

In this section, we perform a brief survey of dynamics of test particles moving within 
potential wells formed along both the equatorial and halo circular orbits. In order to 
do so, we apply several complementary methods of investigation of nonlinear dynamic 
systems. First of all, we construct Poincare surfaces of section which give an overall 
perspective of the phase space dynamics on a given energy hypersurface (for a given 
values of system parameters). For the inspection of individual trajectories, however, we 



prefer to analyse their recurrence plots 34 , which proved to be very useful method in 



our previous work 12 . Besides other properties of recurrence plots (RPs), we highlight 
their ability to clearly distinguish between chaotic and regular dynamics on a short time 
scale, thus reducing the integration time needed for the analysis. We also investigate 
intrinsic frequencies of the orbits employing the rotation number [35] that allows us to 
detect and locate resonances of the system. Recently, the methods of Poincare surfaces 
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Figure 8. Positions of circular halo photon orbits in Bonnor spacetimes. Orbits can 
occur only at radii r > 3.333 and in spacetimes with \b\ > 1.361. 



of section and recurrence plots were very well applied, e.g., in the papers 36 and 37 . 




Figure 9. Off-equatorial trajectories of charged test particle with L = —2.356 a and 
q = 5.581 bounded in halo lobes of the Bonnor spacetime with b = 1 a. In the left panel 
we present a stereometric projection of two trajectories: the upper one with E — 0.8169 
shows ordered motion while with the higher energy the dynamics acquires properties 
of deterministic chaos (bottom trajectory with E = 0.8182). We plot the poloidal 
(r, 6) projection of these trajectories along with several iso-contours of the effective 
potential in the right panel. Both particles were launched at r(0) = 6, 6(0) = ir/3 with 
u r (0) — 0. Grey color marks r — surface in both plots. 

We are dealing with the autonomous Hamiltonian system ^ describing the motion 
of electrically charged test particle around a massive center with dipole magnetic field 
described by the Bonnor metric ([T]). Due to the axial symmetry and stationarity of 
such a setup, we have two obvious constants of motion E and L, related to cyclic 
coordinates t and 0, respectively. Besides that, the value of super-Hamiltonian (|8]) is 
also conserved, which corresponds to the conservation of the particle rest mass. The 



Regular and chaotic orbits near a massive magnetic dipole 



15 



motion of a given particle thus generally occurs on the 3-dimensional hypersurface in the 
phase space. If there is, in addition, another independent integral of motion, analogical 
to the Carter's fourth constant in the Kerr-Newman spacetime (4,38 , the system is 
completely integrable. In such a system, we could find no traces of the deterministic 
chaos at all and, each phase space trajectory would be bound to the 2-dimensio nal 
surface of 3-torus. On the other hand, if there is no extra constant of motion, the 
system is nonintegrable, which means that regular and chaotic orbits coexist in the 
phase space. In general, the amount of chaos depends on the strength of nonintegrable 
perturbation of the originally integrable system. In our case, however, such distinction 
between integrable and nonintegrable parts of the Hamiltonian is not obvious, because 
even if we set b = 0, we do not obtain spherically symmetric Schwarzschild spacetime in 
which the motion of test particle is known to be integrable. Moreover, we have two types 
of possible perturbations: the magnetic field parametrized by b, which itself affects the 
motion of uncharged particles, and the electromagnetic force (proportional to charge q) 
exerted on the particle in the general case q ^ 0, b ^ 0. 

We construct Poincare surfaces of section from the intersection points of trajectories 
with the given 9 = 8 sec = const plane. Instead of plotting the points in pair of canonical 
variables (r, ir r ), we translate the momenta into more intuitive quantity u r and, plot the 
(r, u r ) pairs instead, which does not affect the interpretation of Poincare sections. Given 
dimensionality of our system allows us to distinguish between the regular and chaotic 
orbits easily. The regular orbit bound to the surface of 3-dimensional torus in the phase 
space forms a narrow curve (with zero width, in principle), while the chaotic trajectory 
fills a nonzero surface on the section. In the integrable system, one finds regular orbits 
only. Both types of dynamics, however, are generally present in the nonintegrable case. 

Recurrence plots (RPs) are introduced as a tool of visualizing recurrences of a 



trajectory x(t) in the phase space 39 . The construction of a RP is straightforward 
regardless of the dimension of the phase space. First, we evaluate the binary values of 
the recurrence matrix Ry, which can be formally expressed as 

R ii (e) = e(e-||x(i)-xO')ll), i,j = l,...,N, (29) 

where e is a pre-defined threshold parameter, the Heaviside step function, and N 
specifies the sampling frequency. The sampling frequency is applied to the time segment 
of the trajectory x(t) under examination. There is, however, no unique prescription for 
the appropriate definition of the phase space norm | . | and, simple LP norms are often 
applied. Particularly, L 1 , L 2 (Euclidean norm) or L°° (maximum norm) are the typical 
choice. Nevertheless, we employ a slightly different method (see [12] for details on how 
we measure distances in the phase space for the detection of recurrences). The recurrence 
plot is obtained by direct visualizing of the matrix Ry - black dot is assigned in place 
of unity values, white dot reflects zero values of the matrix elements. The interpretation 



of a resulting plot along with its possible quantification is described in details in 34 



For the analysis of underlying frequencies of the orbit, we use the rotation number 
v. Its evaulation requires, at first, to localize the central fixed point u in the Poincare 
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surface of section which is found in the common center of the primary tori. Being 
equipped with the set of section points of a given trajectory, we calculate the angle 
i?j = Z(7j, between the radius vectors (defined with respect to uq) of each pair of 
succeeding section points Ui and Ui + i. The rotation number v is defined as 

1 - 

v= lim — -> (30) 
1=1 

We actually compute the mean of the finite number of values instead of the limit. 

In the following, we plot the rotation number as a function of the initial value of the 
particle radial coordinate, which reflects its position with respect to the central periodic 
orbit Uq. In such a plot we generally expect to find: i) smooth non-constant segments 
corresponding to regular orbits bound to primary tori with irrational frequency ratios 
(quasiperiodic orbits), ii) plateaus of various width reflecting resonant islands of various 
prominence, iii) 'stochastic' segments resulting from chaotic zones, since the rotation 
number is not well defined for the chaotic orbits (the limit does not converge). Secondary 
islands of the stability are predicted by Poincare-Birkhoff theorem to appear under a 
slight perturbation of the integrable system. They form chains of multiplicity given by 
the denominator of the rotation number, which is rational in this case. Each island 
consists of elliptic fixed point nested in the second order tori. The chain is completed 
by the same number of unstable (hyperbolic) fixed points found between each pair of 
islands. Unstable fixed point appears as an inflection occurring at the same level of 
v as the corresponding plateau in the v{f) plot. However, resonant chains of islands 
of a single multiplicity also appear in some integrable systems. On the other hand, 
the presence of at least two types of islands (with different multiplicities) is regarded 
as strong indication of non-integrability, since no analytic integrable system with this 



property is known so far 35 



6. Motion in equatorial lobes 

In order to compare dynamic properties of the system in all three cases (particle in non- 
magnetized 6 = spacetime, uncharged particle in magnetized spacetime and general 
case q ^ 0, 6 ^ 0), we first investigate the motion in equatorial potential wells, since 



there are no circular halo orbits for 6 = 0. In figure 10, we present series of Poincare 
surfaces of section with 6 sec = tt/2 along with the corresponding plots of a rotation 
number as a function of initial value of radial coordinate. 

Upper left panel shows the Poincare surface of the test particles trajectories when 
the magnetic dipole is switched off by setting 6 = 0. We observe perfectly ordered 
motion with no traces of secondary fixed points nested in Birkhoff islands nor the chaotic 
orbits. Such simple pattern on the section is characteristic for integrable systems. The 
integrability conjecture is also supported by the behaviour of rotation number, which is 
smooth and non-constant throughout the lobe. While we show roughly 50 trajectories on 
the section, the v-plot is constructed from around 1000 of them, which allows much more 
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Figure 10. Comparison of dynamics in equatorial lobes (# m i n = 9 sec = 7r/2). Left 
panels show equatorial Poincare surfaces of section while the corresponding rotation 
numbers of trajectories are presented in the right ones. Common parameters of the 
trajectories in the top panels arc E = 0.951, L = 7.2058 a, q = and b = for 
which the equatorial potential minimum appears at r min = 15 a with V m i n = 0.9494. 
Middle panels show Poincare section and rotation numbers for particles with E — 0.94, 
L = 6.1076 a, q = and b = 2.8535 a which brings the stable circular orbit to 
?~min = 10 a with V m in = 0.9234. Bottom panels are plotted for E = 0.818, 
L = -2.7277a, q = 4.7181 and b = 1 a (while V min = 0.8165 at r min = 6a). 
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Figure 11. Dynamics of test particles in non-magnetized Bonnor spacetime (b = 0) 
near the throat of the equatorial lobe (6> m ; n = 9 SCC = it/2). Left panel shows 
equatorial Poincare surface of section, while the corresponding rotation number v(r) 
of trajectories are presented in the right one. Common parameters of the trajectories 
in both the panels are E = 0.9522, L — 7.2058 a, q = and b — for which the 
equatorial potential minimum appears at r m j n = 15 a with V m i n = 0.9494. Unlike 



the case presented in the upper panels of figure 10 here the potential lobe is opened 
allowing the particles to fall onto the horizon. The island of stability observed in the 
surface of section corresponds to the v = 6/7 resonance. 



detailed inspection for the possible presence of tiny chaotic domains or faint resonances. 

If we perturb the system by the magnetic field, however, the chain of Birkhoff islands 
develops. Although we know that there are some integrable system with resonant islands 
of single multiplicity 35 , here its presence arouses suspicion of nonintegrability, since 
none were present for 6 = 0. Indeed, in the following (see figure 12), we observe chaotic 
motion in this setup (q = 0, b ^ 0), which is irrefutable evidence of being nonintegrable. 

Finally, in the bottom panel of figure 10 , we introduce the charge of the test particle. 
We choose such a combination of parameters that leads to the same value of ratio 
^min/^h, as it is acquired in the previous uncharged case. This makes the two cases 
better comparable and the effect of the newly introduced electromagnetic forces more 
distinguishable. In the last surface of section, we really observe much more complex 
patterns compared to those of uncharged particles. KAM curves of quasiperiodic orbits 
are present as well as several Birkhoff chains of islands corresponding to the resonances 
of intrinsic frequencies of the system. These are interwoven with pronounced chaotic 
layers. Such a picture is typical for considerably perturbed system far from integrability. 

The presented sequence of Poincare surfaces of section suggests the following: the 
motion of the test particles in the Bonnor spacetime without magnetic dipole (b = 0) 
appears perfectly regular and therefore it might be integrable, while the magnetic 
parameter b acts as a nonintegrable perturbation here (indicated by the chaos observed 
in figure 12). With the introduction of ionised test particles, however, the system departs 
considerably farther away from the integrability. 
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The conjecture of integrability of test particle motion in the vicinity of non- 
magnetized Bonnor's solution (which can actually be identified with 5 = 2 member 



of the class of exact Zipoy-Voorhees solutions 40,41]) deserves further examination. 
We first check throughout the parameter space to verify that the dynamics in the closed 
potential lobes always follows the simple regular pattern observed in the upper panels 



of figure [TOj Detailed numerical study revealed no traces of chaos and would therefore 
support the claim of integrability Besides these numerical indications, we also recall 
that Bonnor obtained the static magnetic dipole solution fl8] by transforming the Kerr 



stationary metric in which the test particle motion is known to be integrable 38 . 

However, if we turn our attention from the closed potential lobes to the vicinity 
of saddle points (around which the throat allowing the particle to escape the well and 
fall onto the horizon develops), we eventually observe tiny chaotic layers as illustrated 



in figure 11 These chaotic orbits correspond to those particles which actually leave 
the potential well after certain amount of time. Besides the small fraction of escaping 
particles which exhibit chaotic behavior, the majority of orbits in the lobe remains stable 
and regular. Nevertheless, the presence of chaotic orbits represents a conclusive evidence 
of nonintegrability of a given system. Coincidentally, the nonintegrability of motion in a 



general Zipoy-Voorhees spacetime was very recently shown by Lukes-Gerakopoulos 42 
where the issue was treated in detail. Therefore, we end the discussion here and 
summarize that although the dynamics in the non-magnetized Bonnor spacetime (i.e. 
Zipoy-Voorhees with 5 = 2) is typically regular (we actually found no chaotic orbit in 
closed lobes), the underlying system is not integrable. 

7. Motion in halo lobes 

In the following, we compare dynamics in off-equatorial potential wells for uncharged 



(figure 12) and charged particles (figure 13). For the sake of better comparability, these 



are both chosen to have ^min = 9 sec = tt/3 and equal value of r min /rh. Both series show 
Poincare surfaces of section of particles being launched from the vicinity of off-equatorial 
potential minima (in which the circular halo orbit resides) differing in energy E, which 
governs the size and shape of the lobe (E sets the level at which the effective potential 
surface is being intersected). 

First panels of both figures show the situation if we set the energy E to correspond 
to a small halo lobe, i.e., the energy of the particles is only slightly higher than that 
of corresponding circular halo orbit. In such a case, both the uncharged and charged 
particles show ordered motion - Poincare sections are filled mainly with KAM tori and 
only one deformed resonant island appears. No traces of chaotic dynamics are present 
at all. 

When we increase the energy to the level corresponding to a large halo lobe which 
almost touches its symmetric twin below the equatorial plane, we obtain considerably 



more complex picture on the surfaces of section (upper right panels of figures 12 and 



13). In the uncharged case, we observe well pronounced Birkhoff chains of resonant 
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Figure 12. Poincare surfaces of section (6 SCC = 7r/3) for uncharged particles 
moving within halo lobes in the Bonnor spacetime with b = 5.9771a. Particles with 
L = 3.6743 a are launched from a vicinity of the off-equatorial potential well minimum 
(7"min = 10 a, 9 mm — 9 SCC — 7r/3 and V m - ln — 0.8717) with various values of energy. The 
upper left panel shows the section for the level E — 0.8718 (small halo lobe), in the 
upper right we set E = 0.873 (large halo lobe), E = 0.8739 produces cross-equatorial 
lobe which just emerged from symmetric halo lobes (bottom left panel), while with 
E = 0.88 we obtain the large cross-equatorial lobe which almost opens. 



islands to dominate the scene over receding KAM tori. On the other hand, considerable 
fraction of charged particles under analogous circumstances shows chaotic behaviour. 
Islands of stability are still present, but original KAM tori already disintegrated. 

Once the energy is further increased, the two symmetric halo lobes merge into 
one cross-equatorial region. The dynamics of uncharged particles changes dramatically, 
KAM tori mostly disappear and, only several islands of stability survive the onset of 
chaos (bottom left panel of figure 12 ). However, such abrupt change of the size and shape 
of potential lobe does not considerably affect the motion of charged particles for which 
chaos already appeared inside the halo lobe. The bottom left panel of figure 13 reveals 
that the only difference is that few more islands of stability decompose. Increasing the 
cross-equatorial lobe to its maximal spatial extent causes decomposition of most islands 
of stability, though some of them resist. Surprisingly, more stable islands are then found 
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Figure 13. Poincare surfaces of section (9 SCC — 7r/3) for charged particles with 
q = 0.1259 moving within halo lobes in the Bonnor spacetime with b = 4.5393 a. 
Particles with L — —3.5486 a are launched from a vicinity of the off-equatorial potential 
well minimum (r min = 8 a, 9 min — 9 SCC — n/3 and V mm = 0.8475) with various values 
of energy. The upper left panel shows the section for the level E — 0.8477 (small halo 
lobe), in the upper right we set E = 0.8495 (large halo lobe), E = 0.8496 produces 
cross-equatorial lobe which just emerged from symmetric halo lobes (bottom left panel) 
while with E — 0.851 we obtain the large cross-equatorial lobe. 



for charged particles (bottom right panel of figure 13), while for the uncharged particles 
we observe only few traces of stable motion in the chaotic sea (bottom right panel of 
figure 12). 

We conclude that the energy of particle acts as a trigger for chaotic motion, which 
gradually shifts the dynamics from regular to chaotic. In particular, we confirm that 
motion of uncharged particles in Bonnor spacetime is nonintegrable, which was not clear 
enough from the discussion of equatorial lobes (figure 10). More strongly perturbed 
motion of ionised particles appears to be more prone to chaos, as it shows well inside 
the off-equatorial lobe, while uncharged particles do not become largely chaotic before 
merging both symmetric lobes, which allows them to cross the equatorial plane freely. 

To illustrate the continuous transition from ordered to chaotic dynamics, we pick 
a particular trajectory of charged particle and plot series of Poincare sections along 
with corresponding recurrence plots in figure 14 We launch the particle from the locus 
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Figure 14. Poincare surfaces of section (0 SCC = 7r/3) and corresponding recurrence 
plots for charged particle with q = 5.581 moving within halo lobes in the Bonnor 
spacetime with b = 1 a. Particle with L = —2.356 a is launched just from the locus 
of the off-equatorial potential well minimum (r m - ln = 6 a, 9 mul = 9 SCC = n/3 and 
Vmin = 0.81675) with various values of energy (from the upper left to the bottom 
right: E = 0.8168, E = 0.818, E = 0.8182, E = 0.8183, E = 0.819 and E = 0.8198). 
First three couples of plots show the situation in the halo lobe, while bottom plots 
reveal the dynamics after merging of the lobes. A decisive surface of section cannot be 
constructed for the particle in a opened lobe (E = 0.8198) as it escapes after several 
intersections with the surface, while the corresponding RP shows the chaotic nature of 
the motion unambigiously. Unlike the previous figures, we plot all intersection points: 
downward crossings with u e > (black dot) as well as those resulting from upward 
crossings with u 6 < (red dot) on the Poincare surfaces. 



Regular and chaotic orbits near a massive magnetic dipole 



23 



of the off-equatorial potential minimum with various values of energy E, while other 
parameters remain fixed. The sequence begins with the energy corresponding to a small 
halo lobe where we observe ordered motion manifested by narrow curves on the surface 
of section and, simple diagonal pattern of the recurrence plot (upper left panels of figure 



14). Increasing the energy, however, gradually shifts the dynamics towards deterministic 



chaos - trajectory becomes more and more ergodic as it spans larger fraction of given 
energy hypersurface in the phase space. 

Recurrence plots help to reveal abrupt change of dynamics which occurs when 
the both symmetric halo lobes merge in a single cross-equatorial region. Indeed, we 
observe that RP attains typical chaotic pattern of disrupted diagonal lines when the 
lobes merge (compare upper right couple of panels to bottom left couple in figure 14). 
Bottom middle panels corresponding to large cross-equatorial lobe confirm the full onset 
of chaotic dynamics as the trajectory ergodically fills all allowed region on the surface 
of section, and diagonal segments of RP further disintegrate. The bottom right panel 
corresponds to the opened cross-equatorial lobe from which the particle escapes after 
certain amount of time which is too short for a construction of a decisive Poincare surface 
of section, though, it is long enough to create an unambiguous recurrence plot of chaotic 
nature. We note that the integration time, which was necessary for the construction of 
presented surfaces of section, was typically around 100 times longer than that we needed 
for plotting corresponding RPs. 

The role which the quadrupolar deformation (|7|) plays in inducing the chaotic 
behavior has already been studied in various context in number of works. Gueron and 



Letelier first considered the problem in the framework of classical Newtonian gravity 43 



and, then also for exact solutions to Einstein equations [44]. In the latter paper the 
authors study both the static case of 'distorted' Schwarzschild black hole as well as a 
rotating black hole deformed by a quadrupolar term. The solution they employ in the 
study of a static source is a special case of solution [45] that incorporates the whole 
series of higher multipoles. The rotating object is described by a metric obtained as 
a two-soliton solution using the Belinsky-Zakharov method [46]. From the numerical 
study of test particle trajectories in all these systems they conclude that only the prolate 
quadrupole deformation (Q > 0) induces chaotic dynamics. Since they find no traces 
of chaos for oblate deformations (Q < 0), they suggest that the system could become 



integrable in this case 44 



The issue has been revisited by Dubeibe et al. 47 who employ the solutions which 
are in a close relationship to Bonnor's massive dipole which we analyze in this paper. 
Namely, they study a test particle dynamics in Tomimatsu-Sato 5 = 2 spacetime 48 and 
its generalization found by Manko et al. 30 . In the static limit, the former one reduces 
to Zipoy-Voorhees 5 = 2 spacetime which is identical to the non-magnetized Bonnor 
spacetime obtained by setting b = in ([IJ. On the other hand, the solution by Manko 
et al. generalizes the Bonnor's one by considering rotation and quadrupole moment as 
extra free parameters, while Bonnor's solution is static and quadrupole moment is fixed 
by parameters a and b ([7]). The analysis by Dubeibe et al. leads to the conclusion 
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that there are no chaotic orbits in Tomimatsu-Sato 5 = 2 spacetime (whose deformation 
is always oblate) and the system might be fully integrable. On the other hand, the 
Manko et al. solution admits chaos only for oblate deformations. The last point has 
been contradicted by Han 49 who argues that the chaotic regions found by Dubeibe 
et al. for oblate deformations are not astrophysically relevant, since they only appear 
as narrow zones so close to the center that they would actually collide with the star. 
Han's analysis rather shows that motion around oblate Manko et al. star is typically 
regular, while the prolate deformation induces chaotic dynamics in large regions outside 
the star. 

Both Dubeibe et al. and Han considered the Manko et al. solution without the 
magnetic dipole, being interested in a real deformation of the mass rather than additional 
effective quadrupolar deformation caused by the electromagnetic field. Han only briefly 
comments that including the magnetic dipole does not change the global dynamics and, 
in particular, that it does not bring chaos in oblate case. Only neutral test particles 
have probably been considered. 

Our results cannot be directly contrasted with the works mentioned above, since 
we mainly study the dynamics in a slightly different setup (charged test particles and 
nonzero magnetic dipole). Nevertheless, we observe that in the case of neutral particles 
and oblate effective quadrupole moment ([7]), i.e., for < |fe| < a, the dynamics is 
typically regular as we found no chaos in closed lobes under these circumstances. The 
chaos was present only in very narrow zones formed by escaping orbits in opened lobes 
(see figure 11). However, these zones were found well above the horizon of the central 
body. From this we conclude that although the chaos was extremely rare, the system 
is not integrable even for 6 = (see also recent analysis by Lukes- Gerakopoulos [42]). 
On the other hand, if we consider charged particles around Bonnor's source with oblate 
deformation < \b\ < a, we find large regions of chaos easily. Higher values of dipole 
parameter |6| > a leading to prolate deformation (Q > 0) admit chaotic dynamics 
for both neutral and charged test particles. However, the role of effective quadrupole 
moment Q as a trigger for chaos is not clear from our results, since charged particles 
exhibit the chaotic motion even for Q = (|6| = a). 



8. Conclusions 



Our aim with this paper was to study and compare different approaches towards the 
properties of motion near a magnetized compact body, taking into account the interplay 
of gravitational and electromagnetic effects that act on electrically charged matter. To 
this end we adopted the example of Bonnor's solution for the magnetic dipole as a 
spacetime metric that is given in an analytical (suitably simple) form, while preserving 
the essential features of a massive dipole. Our results are described in several sections 
which help us to proceed through different aspects in a systematic manner. 

Firstly (section 3), we investigated properties of equatorial circular motion and, 
determined the regions of radius accessible for different combinations of the spacetime 
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parameters and the particle constants of motion. Then we explored circular motion 
outside the equatorial plane (section 4). This allowed us to reveal the formation of 
the lobes of stable circulation, which typically occurs when the dipole-type magnetic 
field and the central gravitational attraction act simultaneously on charged particles. 
In comparison with the previously studied cases of halo circular orbits circling around 
compact objects with the presence of the dipole-type magnetic fields, the Bonnor's 
massive dipole permits the existence of even circular halo geodesies. 

In the following parts of the paper (sections 5-7), we presented mainly a numerical 
study of the transition between regular and chaotic motion. We analysed the motion 
in the equatorial potential wells (section 6) in three different cases, namely motion in 
non-magnetized spacetime with 6 = 0, motion of uncharged particle on the magnetized 
background (q = 0, b ^ 0) and dynamics in the general case q ^ 0, b ^ 0. Then, we 
again proceeded to the rich variety of properties of motion in the off-equatorial potential 
wells (section 7). 

Our results show that without magnetic field the system hosts mostly regular orbits 
and the dynamics of test particles resembles fully integrable systems. However, further 
numerical inspection revealed that chaotic orbits are also present in this setup proving 
the system nonintegrable. Then we observed that the magnetic parameter b introduces 
a profound perturbation of the dynamics. Moreover, charge of the particle acts as an 
extra perturbation, which shifts magnetized system even farther from the integrability. 
We also studied the role of particle energy E on the degree of chaos found in the 
system, concluding that it acts as a trigger for the chaotic motion. As the energy 
is gradually increased, the system undergoes a continuous transition from the regular 
behaviour to the chaotic dynamics, being almost fully ergodic on the given hypersurface. 
We illustrated such transition by means of Poincare surfaces of section and recurrence 
plots. The role of quadrupole moment has also been discussed. 

Because of the nature of Bonnor spacetimes, the contribution of presented results 
is more or less in theoretical level. However, the motion of charged particles within 



the off-equatorial lobes itself represents also very interesting astrophysical problem 10 



being, e.g., an essential step in modelling of circling off-equatorial structures around 
relevant compact objects. Such a work is being finished these days as well, describing 



electrically charged perfect fluid tori 50 'levitating' above the equatorial plane in 



spherically symmetric and dipolar magnetic fields. 
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